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velocity dependent randomization traffic cellular automaton 

Sven MaerivoelQ and Bart De MooiQ 

Department of Electrical Engineering ESAT-SCD (SISTAf , Katholieke Universttett Leuven 
Kasteelpark Arenberg 10, 3001 Leuven, Belgium 
(Dated: February 8, 2008) 

Within the class of stochastic cellular automata models of traffic flows, we look at the velocity 
dependent randomization variant (VDR-TCA) whose parameters take on a specific set of extreme 
values. These initial conditions lead us to the discovery of the emergence of four distinct phases. 
Studying the transitions between these phases, allows us to establish a rigorous classification based 
on their tempo-spatial behavioral characteristics. As a result from the system's complex dynamics, 
its flow-density relation exhibits a non-concave region in which forward propagating density waves 
are encountered. All four phases furthermore share the common property that moving vehicles can 
never increase their speed once the system has settled into an equilibrium. 
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I. INTRODUCTION 



In the field of traffic flow modeling, microscopic traffic 
simulation has always been regarded as a time consum- 
ing, complex process involving detailed models that de- 
scribe the behavior of individual vehicles. Approximately 
a decade ago, however, new microscopic models were be- 
ing developed, based on the cellular automata program- 
ming paradigm from statistical physics. The main advan- 
tage was an efficient and fast performance when used in 
computer simulations, due to their rather low accuracy 
on a microscopic scale. These so-called traffic cellular 
automata (TCA) are discrete in nature, in the sense that 
time advances with discrete steps and space is coarse- 
grained (allowing for high-speed simulations, especially 
when they are performed on a platform for parallel com- 
putation) 

The seminal work done by Nagel and Schreckenberg 
in the construction of their stochastic traffic cellular au- 
tomaton (i.e., the STCA) 0, led to a global adoptation 
of the TCA modeling scheme. In order to generate traffic 
jams, their model needs some kind of fluctuations, which 
are introduced by means of a noise parameter. But one 
of the artefacts associated with this STCA model, is that 
it gives rise to (many) unstable traffic jams Q- A possi- 
ble approach to achieve stable traffic jams, is to reduce 
the outflow from such jams, which can be accomplished 
by implementing so-called slow-to-start behavior. The 
name is derived from the fact that vehicles exiting jam 
fronts are obliged to wait a small amount of time. One of 
these implementations, is the approach initially followed 
by Takayasu and Takayasu Q . They made the stochastic 
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noise dependent on the distance between two consecutive 
vehicles. Another implementation is based on making 
this noise parameter dependent on the velocity of a vehi- 
cle, leading to the development of the velocity dependent 
randomization (VDR) TCA. This TCA model moreover 
exhibits metastability and hysteresis phenomena 0- 

In this context, our paper addresses the VDR-TCA 
slow-to-start model whose parameters take on a specific 
set of extreme values. Even though this bears little di- 
rect relevance for the understanding of traffic flows, it 
will lead to an induced 'anomalous' behavior and com- 
plex system dynamics, resulting in four distinct emergent 
phases. We study these on the basis of the system's re- 
lation between density and flow (i.e., the fundamental 
diagram) which exhibits a non-concave region where, in 
contrast to the properties of congested traffic flow, den- 
sity waves are propagated forwards. We furthermore in- 
vestigate the system's tempo-spatial evolution in each of 
these four phases, and consider an order parameter that 
allows us to track the transitions between them. 

Although some related studies exist (e.g., [EHUiISEj 
[Tll|). the special behavior discussed in this paper has not 
been reported as such, as most previously done research 
is largely devoted to empirical and analytical discussions 
about these TCA models, operating under normal condi- 
tions. To strengthen our claims, we compare our results 
with the existing literature at the end of this paper. 

In section Inl we briefly describe the concept of a traf- 
fic cellular automaton and the experimental setup used 
when performing the simulations. The selected VDR- 
TCA model is presented in section Hill as well as an ac- 
count of some behavioral characteristics under normal 
operating conditions. The results of our various exper- 
iments with more complex system dynamics are subse- 
quently presented and extensively discussed in section 
IIVI after which the paper concludes with a comparison 
with existing literature in section and a summary in 
section IVll 
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II. EXPERIMENTAL SETUP 

In this section, we introduce the operational character- 
istics of a standard single-lane traffic cellular automaton 
model. To avoid confusion with some of the notations in 
existing literature, we explicitly state our definitions. 



A. Geometrical description 

Let us describe the operation of a single-lane traffic 
cellular automaton as depicted in Fig. ^ We assume N 
vehicles are driving on a circular lattice containing K 
cells, i.e., periodic boundary conditions (each cell can be 
occupied by at most one vehicle at a time). Time and 
space are discretized, with AT = 1 s and AX = 7.5 m, 
leading to a velocity discretization of AV = 27 km/h. 
Furthermore, the velocity Vi of a vehicle i is constrained 
to an integer in the range {0, . . . , fmax}, with wniax typ- 
ically 5 cells/s (corresponding to 135 km/h). 
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FIG. 1: Schematic diagram of the operation of a single-lane 
traffic cellular automaton (TCA); here, the time axis is ori- 
ented downwards, the space axis extends to the right. The 
TCA's configuration is shown for two consecutive time steps t 
and t + 1, during which two vehicles i and j propagate through 
the lattice. Without loss of generality, we denote the number 
of empty cells in front of vehicle i as its space gap Qs^ . 

Each vehicle i has a space headway hg. and a time 
headway ht- , defined as follows: 



of all the vehicles, based on a rule set that reflects the 
car-following behavior. 

Most rule sets of TCA models do not use the space 
headway hs- or the space gap gg. , but are instead based 
on the number of empty cells di in front of a vehicle i. 
Keeping equation in mind, we therefore adopt the 
following convention from this moment on: we assume 
that, without loss of generality, for a vehicle i its length 
Li — 1 cell. This means that when the vehicle is residing 
in a compact jam, its space headway hg- = 1 cell and its 
space gap is consequently gs^ = cells (i.e., equivalent to 
the number of empty cells dt in front of the vehicle). This 
abstraction gives us a rigorous justification to formulate 
the TCA's update rules more intuitively using space gaps. 

B. Performing measurements 

In order to characterize the behavior of a TCA model, 
we perform global measurements on the system's lattice. 
These measurements are expressed as macroscopic quan- 
tities, defining the global density fc, the space mean speed 
Vs, and the flow q as: 



q = kvs, (5) 

with N the number of vehicles in the system and K the 
number of cells in the lattice. The above measurements 
are calculated every time step, and they should be 
averaged over a large measurement period Tgjj^ in order 
to allow the system to settle into an equilibrium. 



hs, = Li+ gs,, (1) 
hu = Pi + 9u ■ (2) 

In these definitions, gsi and gt^ denote the space and 
time gaps respectively; Li is the length of a vehicle and 
Pi is the occupancy time of the vehicle (i.e., the time 
it 'spends' in one cell). Note that in a traffic cellular 
automaton the space headway of a vehicle is always an 
integer number, representing a multiple of the spatial dis- 
cretisation AX in real world measurement units. So in 
a jam, it is taken to be equal to the space the vehicle 
occupies, i.e., hg. = 1 cell. 

Local interactions between individual vehicles in 
a traffic stream are modeled by means of a rule set 
that reflects the car-following behavior of a cellular 
automaton evolving in time and space. In this paper, 
we assume that all vehicles have the same physical 
characteristics (i.e., homogeneous traffic). The system's 
state is changed through synchronous position updates 



Correlation plots of these aggregate quantities lead 
to time-independent graphs conventionally called 'fun- 
damental diagrams' J2q| (e.g., Fig.|2|). An important re- 
mark is that is that these diagrams are actually only valid 
under stationary traffic conditions and for homogeneous 
traffic, which is almost never the case. As such, real-life 
traffic does not 'lie' on this curve because it represents the 
'average behavior' of the vehicles in a traffic stream. And 
although we should more correctly refer to our measure- 
ments as points in a certain phase space (e.g., the {k^q) 
phase space), we will still use the terminology of 'funda- 
mental diagram' in the remainder of this paper when we 
are in fact referring to this phase space. 

Note that we also can calculate these macroscopic 
quantities as local averages (i.e., local in time and space), 
but these lead to the same results as for large systems 
and measurements periods [T^ . 
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The previous global macroscopic measurements (den- 
sity, average speed, and flow) from equations ©, Q, 
and 13), can be related to the microscopic equations 
and ||2J) as follows: 

k oc h^^, (6) 

q oc \ (7) 

with hs and ht the average space and time headway 
respectively. Note that with respect to the time gaps 
and time headways, we will work in the remainder of 
this paper with the median instead of the arithmetic 
mean because the former gives more robust results when 
^ti,gti — *■ +0O for a vehicle i. 

All the fundamental diagrams in this paper, were cal- 
culated using systems of 10'^ cells. The first 10'^ s of each 
simulation were discarded in order to let initial transients 
die out; the system was then updated for Tgjjj^ = 10** s. 

For a deeper insight into the behavior of the space 
mean speed Vs, the average space gap 'g^, and the 
median time gap g^, detailed histograms showing their 
distributions are provided. These are interesting because 
in the existing literature (e.g., [TsL IT^ these dis- 
tributions are only considered at several distinct global 
densities, whereas we show them for all densities. Each 
of our histograms is constructed by varying the global 
density k between 0.0 and 1.0, calculating the average 
speed, the average space gap and the median time gap 
for each simulation run. A simulation run consists of 
5 X 10'' s (with a transient period of 500 s) on systems 
of 300 cells, varying the density in 150 steps. Note that 
a larger size of the system's lattice, has no significant ef- 
fects on the results, except for an increase of the variance. 

All the experiments were carried out with our Java 
software "Traffic Cellular Automata", which can be 
found at http://smtca.dyns.cx jl6|. 

III. VELOCITY DEPENDENT 
RANDOMIZATION 

In this section, the rule set of the VDR-TCA model 
is explained, followed by a an overview of the model's 
tempo-spatial behavior and its related macroscopic quan- 
tities (i.e., the fundamental diagrams and the distribu- 
tions of the speeds and the space and time gaps) under 
normal conditions. 

A. The VDR-TCA's rule set 

As indicated before, we focus our research on the VDR- 
TCA model for the implementation of the car-following 
behavior. The following equations (based on |^) form its 
rule set; the rules are applied consecutively to all vehicles 
in parallel (i.e., synchronous updates): 



Rl: determine stochastic noise 

j v,{t - 1) = =^ p' ^Pq, fo\ 
\v,{t-l)>0 =^ p'^p, 

R2: acceleration and braking 

Vt{t) <— mm{vi(t - 1) + l,5si(^ - l),wmax}, (9) 
R3: randomization 

C^it) <p' =^v^{t) ^max{0,i;.(i)-l}, (10) 

R4: vehicle movement 

x,it)^x,{t~l)+v,{t). (11) 

In the above equations, Vi{t) is the speed of vehicle 
i at time t (i.e., in the current updated configuration of 
the system), Umax is the maximum allowed speed, gs^ 
denotes the space gap of vehicle i and Xi G {1, . . . , K} an 
integer number denoting its position in the lattice. In the 
third rule, equation (fTn|l . £^i{t) £ [0, 1[ denotes a uniform 
random number (specifically drawn for vehicle i at time t) 
and p' is the stochastic noise parameter, dependent on the 
vehicle 's speed {pq is called the slow-to-start probability 
and p the slowdown probability, with Po,p [0, 1]). 

In a nutshell, rule Rl, equation ©, determines the 
correct velocity dependent randomization. Note that we 
only consider the case with two different noise parame- 
ters (i.e., po and p), ignoring the more general case where 
we have a noise parameter for each possible speed (i.e., 
Po, PvYaax)- R-uIg R2, equation 0, states that a 
vehicle tries to increase its speed at each time step, as 
long as it hasn't reached its maximal speed and it has 
enough space headway. It also states that when a ve- 
hicle hasn't enough space headway, it abruptly adapts 
its speed in order to prevent a collision with the leading 
vehicle (the model thus exhibits an infinite braking ca- 
pability). The randomization parameter determined in 
equation ||SJ|, is now used in rule i?3, equation (|10|l . to 
introduce a stochastic component in the system: a ve- 
hicle will randomly slow down with probability p' . The 
last rule i?4, equation Hll|l . isn't actually a 'real' rule; it 
just allows the vehicles to advance in the system. 

B. Normal behavioral characteristics 

Depending on their speed, vehicles are subject to differ- 
ent randomizations: typical metastable behavior results 
when Po 3> p, meaning that stopped vehicles have to wait 
longer before they can continue their journey (i.e, they 
are 'slow-to-start'). This has the effect of a reduced out- 
flow from a jam, so that, in a closed system, this leads to 
an equilibrium and the formation of a compact jam. For 
such a typical situation (e.g., po = 0.5 and p — 0.01), the 
tempo-spatial evolution is depicted in Fig. [2 each vehicle 



is represented by a single colored dot; as time advances, 
vehicles move to the upper right corner, whereas conges- 
tion waves move to the lower right corner, i.e., backwards 
in space (note that the orientation of the axes is differ- 
ent than the one Fig.^. We can see an initially homo- 
geneous traffic pattern (one metastable phase) breaking 
down and kicking the system into a phase separated state 
(two phases: free-flowing traffic and a compact jam). 




Time (580 s) 



FIG. 2: A typical time-space diagram of the VDR-TCA 
model. The shown lattice contains 300 cells (corresponding to 
2.25 km), the visible time horizon is 580 seconds. The slow-to- 
start probability po was set to 0.5, the slowdown probability 
p to 0.01, and the global density equal to the system's critical 
density. We can see a phase separation taking place, forming 
a persistent compact jam surrounded by free- flowing traffic. 
The significant decrease of the density in the regions outside 
the jam results from the jam's reduced outflow. 

Studying the (fc,g) fundamental diagram in Fig. 13 
gives us another view of this phase transition. The 
free-flow branch was obtained by initially distributing 
the vehicles homogeneously over the system's lattice, 
whereas they were placed in a compact superjam (i.e., 
all vehicles are 'bunched up' behind each other) in order 
to obtain the congested branch. We can see a 'capacity 
drop' taking place at the critical density, where traffic 
in its vicinity behaves in a metastable manner. This 
metastability is characterized by the fact that sufficiently 
large disturbances of the fragile equilibrium can cause 
the flow to undergo a sudden decrease, corresponding 
to a first-order phase transition in statistical physics. 
The state of very high flow is then destroyed and the 
system settles into a phase separated state with a large 
mega jam and a free-flow zone The large jam will 

persist as long as the density is not significantly lowered, 
meaning that recovery of traffic from congestion thus 
shows a hysteresis phenomenon [T^ . 

Considering the non-stochastic CA-184 (i.e.. Wol- 
fram's rule 184 with Vmax > 1 0])' can calculate the 
critical density (where the transition between the free- 
flowing and congested phases occurs) as ,20j : 



kc - y- . (12) 

^ + wmax 
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FIG. 3: Typical {k,q) fundamental diagram of the VDR- 
TCA (with the slow-to-start probability po set to 0.5 and the 
slowdown probability p set to 0.01). The typical 'reversed 
A' shape is obtained by going from the free-flowing to the 
congested traffic regime and vice versa, resulting in hysteretic 
behavior (i.e., the directions of the arrows); in the flrst case, 
traffic encounters a capacity drop, whereas in the second case, 
the recovery results in a lower value for the flow. 

Setting vxnax = 5 cells/s (this equals 5 x 7.5 x 3.6 = 
135 km/h), we correspondingly get kc = (l-f5)^^ = 0.16. 
Note that this number is dimensionless because kc 
represents the fraction of the maximum density, 
i.e., the jam density kj — {AX)~^ = 133.3 ve- 
hicles/km. Converting the value of the critical 
density to real world measurement units, we obtain 
kc = 0.16 X 133.3 = 22.2 vehicles/km. The associated 
capacity flow qcap = kc ■ Wmax = 135 x 22.2 = 3000 ve- 
hicles/h. Note that this number is an overestimate for 
our system, as an exact formula expressing the capacity 
flow gcap in function of the maximum speed Wmax, 
the slowdown probability p, and the lattice size K, has 
not yet been determined |2l| . If we perform a linear 
regression on the congested branch of Fig. |21 we get the 
speed of the backward propagating congestion waves, 
which is « —13.5 km/h, close to the empirical value of 
-15 km/h found in real-life traffic 

If we look at the distribution of the vehicles' speeds, 
we get the histogram in Fig. 0] Here we can clearly 
see the distinction between the free-flowing and the 
congested phase: the space mean speed remains constant 
at a high value, then encounters a sharp transition (i.e., 
the capacity drop), resulting in a steady declination as 
the global density increases. Note that as the critical 
density is encountered, the standard deviation jumps 
steeply; this means that vehicles' speeds fluctuate wildly 
at the transition point (because they are entering and 
exiting the congestion waves). Once the compact jam 
is formed, the dominating speed quickly becomes zero 
(because vehicles are standing still inside the jam). 
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Global density k Histogram class [speed v in cells/s] 

FIG. 4: The distribution of the vehicles' speeds v, as a function of the global density k in the VDR-TCA (with po = 0.5 and 
p = 0.01). In the contourplot to the left, the thick solid line denotes the space mean speed, whereas the thin solid line shows its 
standard deviation. The grey regions denote the probability densities. The histograms (A) and (B) to the right, show two cross 
sections made in the left contourplot at = 0.1325 and k = 0.4000 respectively: for example, in (A), the high concentration of 
probability mass at the histogram class v = 5 cells/s corresponds to the dark region in the upper left corner of the contourplot. 



Although most of the weight is attributed to this 
zero-speed, there is a non-neghgible maximum speed 
present for intermediate densities. If the global density is 
increased further towards the jam density, this maximum 
speed disappears and the system settles into a state in 
which all vehicles either have speed zero or one (i.e., 
systemwide stop-and-go trafhc). 

Considering the distribution of the vehicles' space 
gaps, we get Fig.|5l because of their tight coupling in the 
VDR-TCA's rule set, the courses of both the space mean 
speed and the average space gaps are similar. Although 
the space gaps are rather large for low densities, at the 
critical density kc they leave a small cluster around an 
optimal value. This value can be calculated using equa- 
tions and ifT^ (noting that L = L ~ 1 because 
we are dealing with homogeneous traffic): 



g^^h,-L^-~L, (13) 

resulting in [g,] = [0.1669-1 - 11 = [4-99] = 5 cells. 
This corresponds to the necessary minimal space gap 
in order to travel at the maximum speed, avoiding a 
collision with the leader. An important observation is 
that no recorded space gaps exist between this cluster 
of five cells and the space gaps of zero cells inside the 
compact jam. This means that there is a distinct phase 
separation taking place once beyond the critical density: 
vehicles are either completely in the free-flowing regime, 
or they are in the compact jam. 

Figure IHl shows the distribution of the vehicles' time 
gaps: when traffic is in the free-flowing regime, time gaps 




Global density k 

FIG. 5: The distribution of the vehicles' space gaps Qs, as 
a function of the global density k in the VDR-TCA (with 
po ~ 0.5 and p = 0.01). The thick solid line denotes the 
average of all the vehicles' space gaps, whereas the thin solid 
line shows its standard deviation. The grey regions denote 
the probability densities. 



are high, but finite. As the critical density is approached, 
the median time gap flrst decreases (because the vehicles' 
speeds remain the same but their space gaps decrease). 
Once beyond the critical density, it increases towards 
infinity because vehicles come to a full stop inside the 
compact jam. Just as with the space gaps, we can also 
observe a small cluster around an optimal value. This 
optimal time gap, is the time needed to travel the dis- 
tance formed by the optimal space gap, at the maximum 
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speed. This means that = Wmax = 1 s. Because 
the previous train of thought is only exactly valid for the 
non-stochastic CA-184 and the VDR-TCA incorporates 
stochastic noise, the median optimal time gap lies some- 
what above the previously calculated value {g^ « 1.2 s). 



use of an order parameter that is able to track the phase 
transitions between the different traffic regimes. 
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Global density k 

FIG. 6: The distribution of the vehicles' time gaps gt, as 
a function of the global density k in the VDR-TCA (with 
po ~ 0.5 and p = 0.01). The thick solid line denotes the 
median of all the vehicles' time gaps. The grey regions denote 
the probability densities. 



IV. MORE COMPLEX SYSTEM DYNAMICS 

Most of the previous research dealt with the study 
of the behavioral characteristics of the VDR-TCA 
model operating under normal conditions. We now 
turn our attention to the specific case in which the 
model's parameters take on extreme values pq <^ p, more 
specifically considering the limiting case where po = 0.0 
and p — 1.0. 

We will first look at the change in tempo-spatial be- 
havior when p is increased towards 1.0, at which point a 
peculiar behavior is established in the system. We study 
the qualitative effects that the VDR-TCA's rule set has 
on individual vehicles, and discuss shortly the prevail- 
ing initial conditions. This is followed by a quantitative 
analysis using the {k,q) and (fc,?7s) fundamental diagrams. 
The former diagram leads us to the discovery of four dis- 
tinct phases, having a non-concave region with forward 
propagating density waves. More elaborate explanations 
are given based on the histograms of the space mean 
speed, the average space gap, and the median time gap. 
This is followed by a detailed analysis of the observa- 
tions of the tempo-spatial behavior of the system in each 
of the four different phases (i.e., traffic regimes). The 
section concludes with a short discussion on the effects 
of different maximal speeds, after which we adopt the 



A. Increasing the stochastic noise p 

Let us first consider the case in which pq = 0.0 and 
where we vary p between 0.0 and 1.0. Figure H shows 
a time-space diagram where we simulated a system con- 
sisting of 300 cells. As time advances (over a period of 
580 s), the slowdown probability p is steadily increased 
from 0.0 to 1.0 (the global density k was set to 0.1667 
which is slightly below the critical density for a system 
with stochastic noise). 
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Time (580 s) with p e [0 ^ 1] 

FIG. 7: A time-space diagram showing a system with a lat- 
tice of 300 cells (corresponding to 2.25 km); the visible time 
horizon is 580 seconds. As time advances, the slowdown prob- 
ability p varies between 0.0 and 1.0 (the slow-to-start proba- 
bility Po was fixed at 0.0). The system's global density k is 
0.1667. 



When p is very low, vehicles can keep driving in the 
free-fiowing regime. As p is increased, small unstable 
jams occur. Increasing p even further, leads to even more 
pronounced jams. An important observation is that the 
propagation speed of these backward moving congestion 
waves increases; whenp reaches approximately 0.5, their 
speed equals 0, meaning that jams stay fixed at a posi- 
tion. Note that as these jams are unstable, they can be 
created or dissolved at arbitrary locations. Finally, as p 
tends to 1.0, we can see the emergence of forward prop- 
agating congestion waves; 'density' is now being carried 
in the direction of the traffic flow. These forward propa- 
gating waves form 'moving blockades' that trap vehicles, 
which in turn leads to tightly packed clusters of vehicles 
that move steadily (but at a much slower pace than in the 
free-fiowing regime). Note that this 'clustering' behavior 
is different from platooning, which typically occurs when 
vehicles are driving close to each other at relatively high 
speeds 22] . In our case, the clusters of vehicles advance 
more slowly. 
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B. Qualitative effects of the rule set 

From now on, wc only consider the limiting case where 
Po — 0.0 and p = 1.0. Let us now study the influence of 
the VDR-TCA's rule set on an individual vehicle i. As- 
suming Vi G {0, . . . , Wmax}, the rules described in section 
1111 Al lcad to the following four general cases: 

• Case (1) with Vi{t — 1) = Q and gsi{t — 1) = 

In this case, the vehicle is residing inside a jam 
and rule R2 (acceleration and braking) plays a 
dominant part: the vehicle's speed Vi{t) remains 0. 

• Case (2) with - 1) = and gs^(t - 1) > 

This situation may arise when, for example, a 
vehicle is at a jam's front. According to rule 
Rl, the stochastic noise parameter p' becomes 
0.0. Rule R2 then results in an updated speed 
Vi{t) — inm{l,gs-{t — 1)} — 1. Because there is no 
stochastic noise present in this case, rule iZ3 does 
not apply and the vehicle always advances one cell. 



From the table it follows that all vehicles 
traveling at Vi < Wmax can neither accelerate 
nor decelerate: the vehicles' current speed is 
kept. 

- Case (4b) with v,{t - 1) > gs^{t - I) 

Just as in the previous case (4a), the stochas- 
tic noise becomes 1.0. Rule R2 now changes 
to Vi{t) ^ gs.(t—l). Because rule i?3 is always 
applied, this results in vi (t) actually becoming 
gs-{t — l) — l instead of just gs-{t — l). So a ve- 
hicle always slows down too much (as opposed 
to solely avoiding a collision with its leader). 

Conclusion: considering the previously discussed four 
general cases, the most striking feature is that, according 
to case (4a), in a VDR-TCA model with po = 0.0 and 
p — 1.0, a moving vehicle can never increase its speed, 
i.e., Vi{t) < v.i{t - 1). 

C. Effects of the initial conditions 



• Case (3) with v,{t - 1) > and gs,{t - 1) = 

Because there is no space in front of the vehicle, 
it has to brake in order to avoid a collision. 
Rule R2 consequently abruptly decreases the ve- 
hicle's speed Vi{t) to 0, as the vehicle needs to stop. 

• Case (4) with Vi{t -l)> Q and gs^{t - 1) > 

This case deserves special attention, as there are 
now two discriminating possibilities: 

- Case (4a) with Vi{t ~ I) < gs^it - I) 

According to rule Rl, the stochastic noise p' 
becomes 1.0. Because the vehicle's speed is 
strictly less than its space gap, rule R2 be- 
comes Vi{t) ^ min{ui(t - 1) -l- l,i'max}- Fi- 
nally, rule iZ3 is applied which always reduces 
the speed calculated in rule R2 (constrained 
to 0) . In order to understand what is happen- 
ing, consider the speeds Vi{t — 1) and Vi{t) in 
the following table: 

Vi{t-1) v^{t) 

vmax — > Umax - 1 

Umax - 1 — > Umax - 1 

Umax - 2 — > umax - 2 

2 — y 2 
1 — y 1 

We can clearly see that the maximum speed a 
vehicle can travel at, is constrained by Wmax 
- 1, which corresponds to fmax — p' uM- 



As already mentioned, in this paper we study the 
limiting case where po = 0.0 and p = 1.0. This case is 
special, in the sense that the system's behavior for light 
densities is extremely dependent on the initial conditions. 

If we calculate the flow q associated with each global 
density k in the range [0, -j], we obtain Fig. |S1 Here we 
show three curves for Vmax S {3, 4, 5} where the system 
was initialized by distributing the vehicles homoge- 
neously over the system's lattice. A fourth curve is also 
shown, corresponding to wmax = 5 where, as opposed 
to the previous initialization scheme, the system was 
initialized by distributing the vehicles randomly over the 
system's lattice. 

After the system has settled into an equilibrium, the 
resulting flow is limited by the fact that the maximum 
speed of a vehicle in the system is always equal to 
Vmax — 1. As proven in section FlVBI vehicles can never 
accelerate, which means that the slowest car in the sys- 
tem determines the maximum possible flow. Therefore, 
if the system (with Wmax = 5 cells/s) is initialized with a 
homogeneous distribution of the vehicles, then all of them 
will travel at Wmax — 1 = 4 cells/s (i.e., the upper branch 
in Fig.|SJ). If however, the system is initialized randomly, 
this has the effect that some vehicles are spaced more 
closely to each other. As a direct consequence of this, all 
vehicles will now travel at f^ax ~ ^ cells/s with wj^iax 
the speed of the slowest vehicle in the system. In the 
extreme case, a single vehicle traveling at = 1 ccll/s 
will cause all other vehicles to slow down to this speed, 
resulting in the lower branch in Fig. |S1 

Note that the effects of different initial conditions (i.e., 
homogeneous versus random distributions of the vehicles) 
are only relevant for light densities. Because the phase 
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FIG. 8: The VDR-TCA (with po = 0.5 and p = 0.01) sys- 
tem's behavior for hght densities has a strong dependence on 
the initial conditions. The {k,q) fundamental diagrams are 
shown for vmax. G {3, 4, 5} with k G [0, The systems have 
been initialized using either a homogeneous or a random dis- 
tribution of the vehicles. 



transitions in this paper occur at moderate to high densi- 
ties, it doesn't matter what kind of initiaUzation scheme 
is used, so from now on we always assume that all the 
vehicles are distributed randomly. 



D. Quantitative analysis 

Whereas the previous sections dealt with the effects 
of the VDR-TCA's changed rule set and the role of the 
initial conditions, this section considers the effects on the 
and {k,Vs) fundamental diagrams (i.e., three phase 
transitions and a non-concave region), as well as on the 
histograms of the space mean speed, the average space 
gap and the median time gap. 



1. Effects on the fundamental diagrams 



ignore the hysteretic behavior of the fundamental dia- 
grams, as it gives no relevant contribution to the results 
of our experiments. 
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FIG. 9: The {k,q) fundamental diagrams of the VDR-TCA 
with a fixed slow-to-start probability po = 0.0 (the slowdown 
probability p is increased from 0.0 to 1.0). In the limiting case 
where p = 1.0, four distinct density regions {I)-{IV) appear. 

As the slowdown probability p is increased from 0.0 
to 1.0, the critical density - at which the transition from 
the free-flowing regime occurs - is shifted to lower values 
(note that the magnitude of the capacity drop also 
diminishes). For a global density k — ^ (i.e., the critical 
density of the non-stochastic CA-184), we can see that 
the speed of the backward propagating congestion waves 
increases, just as was visible in the time-space diagram 
of Fig.[7| The speed of these characteristics is defined as 
dq/dk (i.e., the tangent to the fundamental diagram), 
and when p « 0.5, the sign of the speed is reversed, 
leading to the earlier mentioned forward propagating 
density waves. Furthermore, as is apparent in Fig.|3 for 
high densities there exists a region {IV) in which the 
flow q is only dependent on k and not on the stochastic 
noise p: all the measurements coincide on this heavily 
congested branch. 



Most existing {k,q) fundamental diagrams related to 
traffic flow models, show a concave course (althou gh s ome 
(pedagogic) counter-examples such as 0, S H 1^ llOt ITlj 
exist). Non-concavity (i.e., convexity) of a function / is 
defined as 'ixi,X2 e dom/ | f{^{xi + X2)) < ^/((xi) -t- 
f{x2))- The property of concavity also holds true for 
the VDR-TCA model operating under normal conditions 
(cfr. Fig. 13). However, whenpo ^ P, this is no longer the 
case: in the limit where po = 0.0 and p = 1.0, the fun- 
damental diagram exhibits two distinct sharp peaks. Be- 
tween these peaks, there exists a region where 
the fundamental diagram has a convex shape, as can be 
seen from Fig. |51 In region (///), traffic flow has a ten- 
dency to improve with increasing density. Note that we 



Increasing the slowdown probability p has also an ef- 
fect on the space mean speed in the free-flowing regime. 
As already stated in section ITlI Bl this speed is equal to 
Wmax —p ■ The effect is more pronounced if we look at 
the {k,Vs) fundamental diagram in Fig. EH Here we can 
see that the average maximum speed in the free-flowing 
regime shifts downwards, reaching a value of 4 cells/s 
when p = 1.0. From then on, there are two regions 
(/) and {III) where the average speed remains constant. 
Note that, to be precise, region (/) actually contains a 
small capacity drop at a very low density, but we ignore 
this effect, thus treating region (/) in an overall man- 
ner. Finally, note that the {k,Vs) fundamental diagram 
remains a decreasing function as p — > 1.0. 
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FIG. 10: The {k,Vs) fundamental diagrams of the VDR-TCA 
with a fixed slow-to-start probability po = 0.0 (the slowdown 
probability p is increased from 0.0 to 1.0). 



2. Effects on the histograms 




0.3 0.4 0.5 0.6 0.7 

Global density k 

FIG. 11: The distribution of the vehicles' speeds v, as a, 
function of the global density k in the VDR-TCA (with po = 
0.0 and p = 1.0). The thick solid line denotes the space 
mean speed, whereas the thin solid line shows its standard 
deviation. The grey regions denote the probability densities. 



The distribution of the vehicles' speeds is shown in 
Fig. ^] just as in the {k,Vs) fundamental diagram in 
Fig- El we can clearly observe two 'probability plateaus'. 
In the first region (/), vehicles' speeds are highly concen- 
trated in a small region around Vs — 1 cell/s. As the 
global density increases in region (//), the space mean 
speed declines until it reaches region (///) where the sec- 
ond plateau is met at Vs = 0.5 cell/s. From then on, 
it steadily decreases, reaching zero at the jam density 
{kj = 1.0). Note that in region (/), the standard devia- 
tion is zero, whereas it is non-zero but constant in region 
(///). This means that vehicles in the former region all 
drive at the same speed; in the latter region they drive at 
speeds alternating between and 1 cell/s (i.e., stop-and- 
go traffic). 

Considering the steep descending curve at the begin- 
ning of region (/), we state that although vehicles can 
drive at w = Vmax. — 1 = 4 cell/s under free-flowing 
conditions, they nonetheless all slow down as soon as 
at least one vehicle has a too small space gap. In other 
words, when g^. < wniax, case (4b) from section IIV Bl 
applies and the vehicle slows down. This leads to a chain 
reaction of vehicles slowing down, because vehicles can 
never accelerate, as was pointed out in section llV CI 

Under normal operating conditions (i.e., po S> p), a 
vehicle's average speed and average space gap show a 
high correlation in the congested density region beyond 
the critical density. This can be seen from the similar- 
ity between the histogram curves in Figs. 0] and [5| In 
high contrast with this, is the distribution of the vehi- 
cles' space gaps as in Fig. ^1 which shows a different 
scenario. More or less similar to the vehicles' speeds, we 



can observe the formation of three plateaus of constant 
space gaps in certain density regions. These plateaus are 
located at 2, 1 and cells for density regions (/), (///), 
and {IV) respectively. Note that the average space gaps 
themselves are not constant in these regions, as opposed 
to the space mean speed. This is due to the fact that ve- 
hicles encounter waves of stop-and-go traffic, whereby the 
frequency of these waves increases as the global density 
is augmented. 

Another observation that we can make, is that the 
standard deviation goes to zero at the transition point 
between density regions (/) and (//). This means that, 
as expected, the traffic flow at this point consists of 
completely homogeneous traffic, in which all the vehicles 
drive with the same space gap gs = 2 cells (and as 
already mentioned, with the same speed v = I cell/s). 

Just as with the previous histograms, there also ap- 
pear to be plateaus of concentrated probability mass in 
the distribution of the vehicles' time gaps in Fig. El As 
opposed to the standard behavior in Fig. |H| the concen- 
tration in the first region (/) is more elongated and more 
or less completely flat. This is expected because the ma- 
jority of the space mean speeds and the average space 
gaps remain constant in this region. Furthermore, once 
the first phase transition occurs, the median time gap 

— > -foo as the space mean speeds and average space 
gaps tend to zero. This is expressed as the existence of 
a non-neglibible cluster of probability mass at the top 
of the histogram in Fig. El the concentration is formed 
by vehicles that are encountering the earlier mentioned 
stop-and-go waves, resulting in the fact that their time 
gaps periodically become infinity. 
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FIG. 12: The distribution of the vehicles' space gaps Qs, as 
a function of the global density k in the VDR-TCA (with 
po = 0.0 and p = 1.0). The thick solid line denotes the 
average of all the vehicles' space gaps, whereas the thin solid 
line shows its standard deviation. The grey regions denote 
the probability densities. 
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As the global density k is increased, the system under- 
goes three consecutive phase transitions (between these 
four traffic regimes). The densities at which these transi- 
tions occur, are not just inflection points (such as the 
finite-size effect in the TCA model of Takayasu and 
Takayasu 0), but they signal a new behavior in the sys- 
tem. The next four paragraphs give a detailed account of 
the effects that appear in each density region (i.e., traf- 
fic regime), as well as the transitions that occur between 
them. 



1. Region (1) - free-flowing traffic [FFT] 

For moderately low densities, typical time-space dia- 
grams look as the ones in Fig. E| According to the his- 
togram in Fig. ^2 we can see that the speed of all the 
vehicles is the same (namely 1 cell/s) in this density re- 
gion. Although the standard deviation of the space mean 
speed is zero, this is not the case for the average space 
gap, explaining its rather 'nervous' behavior in Fig. 1121 

As the global density increases, the transition point 
between regions (/) and (//) is reached. At this point, 
each vehicle i has a speed Vi = 1 cell/s, and a space gap 
gs^ = 2 cells. Because Vi < gs^, case (4a) from section 
II V lj| applies. Using equations ^ and ©, we can calcu- 
late the corresponding density: 



(14) 



Although the maximum speed any vehicle can (and 
will) travel at is limited to 1 cell/s, we still call this state 
'free-flowing traffic ' (FFT) because no congestion waves 
are present in the system and none of the vehicles has to 
stop. 



2. Region (11) - dilutely congested traffic [DCT] 



FIG. 13: The distribution of the vehicles' time gaps gt, as 
a function of the global density k in the VDR-TCA (with 
Po = 0.0 and p = 1.0). The thick solid line denotes the 
median of all the vehicles' time gaps. The grey regions denote 
the probability densities. 



E. Typical tempo-spatial behavior 

Studying the {k,q) fundamental diagrams (Figs. 1^1 
and I10|) in the previous section, we saw the emergence 
of four distinct density regions {I)-{IV) as p 1.0 
(formed by the A-symbols). In this section, we discuss 
the tempo-spatial properties that are intrinsic to these 
regions, relating them to the previously discussed 
histograms. 



If we increase the density to k ~ ^ + ^ (i.e., adding 
one vehicle at the transition point), a new traffic regime 
is entered. In this regime, each extra vehicle leads to a 
backward propagating mini-jam of 3 stop-and-go cycles 
(see Fig. I15II . bringing us to the description of 'dilutely 
congested traffic' (DCT). 

In order to calculate the second transition point, we 
again observe the histogram of Fig.^] A vehicle's space 
gap now alternates between 1 and 2 cells in density region 
(//). And because its speed alternates between and 
1 cell/s respectively, the vehicles' motions are controlled 
by cases (2) and (4a) from section lTV Bl This means that 
stopped vehicles accelerate again in the next time step, 
after which they have to stop again, and so indefinitely 
repeating this cycle of stop-and-go behavior. 




FIG. 14: Time-space diagrams of the VDR-TCA model with 
po = 0.0 and p = 1.0. The shown lattices each contain 300 
cells (corresponding to 2.25 km), the visible time horizon is 
always 580 seconds. The upper diagram and lower diagrams 
have densities k = 0.2 and k — ^ respectively, correspond- 
ing to observations in density region (/), labeled free-flowtng 
traffic (FFT). 
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FIG. 15: Time-space diagrams of the VDR-TCA model with 
Po = 0.0 and p = 1.0. The shown lattices each contain 300 
cells (corresponding to 2.25 km), the visible time horizon is 
always 580 seconds. The upper diagram and lower diagrams 
have densities k = ^ + and k = 0.4 respectively, corre- 
sponding to observations in density region (//), labeled di- 
lutely congested traffic (DOT). 



Consider now a pair of adjacent driving vehicles i and 
j; it then follows from equation ^ that hg- = 1 + 1 = 
2 cells and hg^ =1-1-2 = 3 cells. So each pair of vehicles 
'occupies' 5 cells in the lattice, or 2.5 cells on average per 
vehicle. This leads to the second transition point being 
located at: 



=hg =j^ = 0.4. (15) 

As the density is increased towards k(ii)^i^iii), the 
space mean speed decreases non-linearly. At the tran- 
sition point itself, Vs — 0.5 cells/s and the system is now 
completely dominated by backward propagating dilute 
jams. 

3. Region (IIIJ - densely advancing traffic [DAT] 

Adding one more vehicle at the transition point be- 
tween regions (//) and {III), leads to surprising behav- 
ior: a forward moving jam emerges, traveling at a speed 
of 0.5 cells/s (see Fig. I16|l . Another artefact is that the 
cycle of alternating space gaps of 1 and 2 cells, is bro- 
ken with the introduction of a zero space gap at the lo- 
cation of this new 'jam'. When the density reaches the 



third transition point, the system is completely filled with 
these forward moving 'jams' of dense traffic, leading to 
the description of 'densely advancing traffic ' (DAT) . Be- 
cause the available space is more optimally used by the 
vehicles, an increase of the density thus has a (tempo- 
rary) beneficial effect on the global flow measured in the 
system. This kind of behavior of forward moving den- 
sity structures can also be observed in some models of 
anticipatory driving [2^ . 

At the transition point itself, all vehicles exhibit the 
same behavior, comparable to the behavior at the sec- 
ond transition point. Traffic is more dense, as can be 
seen in the distribution of the space gaps in Fig. ^] all 
vehicles have alternatingly = and gs = 1 cell (with 
corresponding speeds of and 1 cell/s respectively). One 
pair of adjacent driving vehicles i and j thus 'occupies' 
+hs- = (l-l-O)-l-(l-l-l) = 3 cells, or 1.5 cells on average 
per vehicle. The third transition point thus corresponds 
to: 



=h, = = g- (16) 



I 




FIG. 16: Time-space diagrams of the VDR-TCA model with 
po = 0.0 and p = 1.0. The shown lattices each contain 300 
cells (corresponding to 2.25 km), the visible time horizon is 
always 580 seconds. The upper diagram and lower diagrams 
have densities k — 0.4+ i and fc = | respectively, correspond- 
ing to observations in density region {III), labeled densely 
advancing traffic (DAT). 
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FIG. 17: Time-space diagrams of the VDR-TCA model with 
Po ~ 0.0 and p = 1.0. The shown lattices each contain 300 
cells (corresponding to 2.25 km), the visible time horizon is 
always 580 seconds. The upper diagram and lower diagrams 
have densities k = ^ + ^ and k = 0.85 respectively, corre- 
sponding to observations in density region {IV), labeled heav- 
ily congested traffic (HCT). 



4-. Region (IV) - heavily congested traffic [HCT] 

Finally, as the system's global density is pushed to- 
wards the jam density, each extra vehicle introduces at 
any point in time a backward propagating jam, consist- 
ing of a block of 5 consecutively stopped vehicles (see 
Fig. El • The pattern of stable stop-and-go traffic gets 
destroyed, as vehicles remain stopped inside jams for 
longer time periods (cfr. density region {IV) in the dis- 
tribution of the time gaps in Fig. lT^ . leading to the de- 
scription of this traffic regime as 'heavily congested traffic ' 
(HCT). 

F. Varying the maximum speed 

Up till now, the paper focused on VDR-TCA 
model's behavior when po = 0.0 and p = 1.0, for 
Vmax — 5 cells/s. We can now ask the question: "What 
would be the effect if we consider a lower maximum 
speed ?" 

Considering the {k,q) fundamental diagrams in Fig. 1181 
we can see that a lower maximum speed has the signifi- 
cant effect of seemingly losing the first two transitions 
points -♦(//) and (///). Put more correctly. 



when Wmax = 1 cell/s, the ffist transition point fc(/)^(/7) 
gets drastically shifted towards a higher density of |. A 
side-effect is that the second and third transitions points 
vanish. This effect is also visible in the {k,Vs) fundamen- 
tal diagrams in Fig.^] Note that we have set p = 0.9 in 
order to make the differences between the fundamental 
diagrams more pronounced. 

A decrease of the maximum speed thus has only a sig- 
nificant effect in the limiting case where Wniax = 1 cell/s, 
leading to a less rich dynamical behavior in the system. 
Any other maximum speed Wmax > 1 cell/s gives rise to 
different traffic regimes. Note that all the branches in 
the congested regions of the fundamental diagrams co- 
incide once the ffist transition between regions (/) and 
(//) occurs at /e(/)^(/7) = |. 

G. Tracking the phase transitions 

In order to track the transitions that occur between 
the different regimes {I)-{IV), we look for a suitable or- 
der parameter that expresses a qualitative behavior that 
changes between these regimes. In the following para- 
graphs, we describe the use of two such parameters: the 
first one is based on correlations between neighboring 
cells, whereas the second one is based on a comparison 
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FIG. 18: The {k,q) fundamental diagrams of the VDR-TCA 
with fixed slow-to-start and slowdown probabilities of po = 
0.0 and a, p = 0.9 respectively. The maximum speed was set 
at 1, 2, and 5 cells/s. 
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FIG. 19: The {k,Vs) fundamental diagrams of the VDR-TCA 
with fixed slow-to-start and slowdown probabilities of po = 
0.0 and a p = 0.9 respectively. The maximum speed was set 
at 1, 2, and 5 cells/s. 



of the difference between the system's global and local 
densities. 



1. Nearest neighbor correlations 

One example of such a suitable measure, is the nearest 
neighbor order parameter Mi , which gives the time-space 
averaged density of those vehicles with speed v = Q cells/s 
which had to brake due to the next vehicle ahead |2(ll24l |. 

If we define the occupancy rji of cell i as: 



rji = 1 4=> cell i is occupied by a vehicle, 
i]i — 4^ cell i is empty, 



(17) 



then we can define the order parameter Mi as the fol- 
lowing function that captures the density correlation be- 
tween pairs of nearest neighbors (at a certain global den- 
sity k): 



Ml = 



(18) 



Sim 



with riK+i = ^Ji- Note that equation H18|l is calculated 
as an average over the total simulation period T'gjj^. 
Figure 1201 shows the results from an experiment for 
calculating Mi when K = 300 cells, Umax = 5 cells/s 
and Tg^j^ — 10^ s (with a transient period of 10^ s). 
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FIG. 20: The Mi nearest neighbor order parameter for the 
two systems with {po,p) = (0.5,0.01) and (po,p) = (0.0,1.0) 
respectively. The system size was 300 cells and the total sim- 



ulation period T„j 



10^ s. 



For the normal case where po = 0.5 and p — 0.01, we 
can see that Mi is clearly able to detect the transition 
at the critical density: the order parameter shows a sud- 
den significant increase. Although the VDR-TCA model 
exhibits metastability, it has only two distinct phases, 
namely free-flowing and jammed. These two phases are 
detected by the order parameter as it is constantly zero 
in the former, whereas it's a linearly increasing function 
of the density in the latter. As stated, in the free-flowing 
regime, the nearest neighbor correlations are always zero 
because all the vehicles have large enough space gaps 
(i.e., at least the cell immediately in front of a vehicle is 
always empty) . 

If we consider the behavior of the order parameter Mi 
in the second case, where po = 0.0 and p = 1.0, we 
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see that it is not able to detect the first transition at 
^(1)^(11) — \- The reason for this can be determined by 
looking at the tempo-spatial evolutions in Figs. 1141 and 
[131 Here we can see that, as described in sections llV E II 
and IIV E 21 each vehicle always has at least one empty 
cell in front of it, leading to a nearest neighbor density 
correlation of zero. 

The order parameter is however succesful in detecting 
the second transition at fc(/7)-+(///) = 0.4. The third 
transition at fc(7/7)_,(7V') = | is detected, but in a 
non-intuitive manner: only the slope of the linear 
dependency between Mi and k changes when going from 
region (///) to region {IV). 

We conclude this paragraph by stating that the Mi 
order parameter is able to capture some transitions be- 
tween different traffic regimes, but not all of them. Fur- 
thermore, it seems that the ability of the parameter to 
distinguish between two regimes is qualitively inferior: 
this can be witnessed by the weak behavioral change at 
the transition point fc(///)^(/y). 



2. A measure of inhomogeneity 

A more suitable order parameter that can accurately 
track the transitions between the different regimes, is 
the order parameter M2 defined by Jost and Nagel |25| . 
They partition the road into L equal segments of integral 
length K/L. After calculating the local densities fc/. for 
all these segments i, they consider the variance of the 
local densities as their order parameter : 



Mi 



(19) 



Sim 



In equation (|19|l . E[ki\ denotes the expected value of 
the local density, which is equal to the global density k 
because we consider a closed system. The local densities 
are calculated as: 



Tgjjj^. Care should be taken not to prematurely average 
the local densities over more than one time step, as this 
would significantly deteriorate the tracking quality of 
the order parameter: the longer the aggregation period 
for calculating the local densities, the more each segment 
'sees' the densities recorded in the other segments as 
congestion waves are traveling around the closed system. 



In Fig.EIlwe can see the results for different slowdown 
probabilities p (the slow-to-start probability po was fixed 
at 0.0). The system size K was set to 300 cells, Wmax — 
5 cells/s and Tgjj^ = 10^ s (with a transient period of 
10^ s). 
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FIG. 21: The M2 order parameter which measures the vari- 
ance of the local densities. In each diagram, the slow-to- 
start probability po was fixed at 0.0, whereas the slowdown 
probability p was varied between 0.1 and 1.0 (the system size 
was 300 cells and the total simulation period Tgiiu = 10^ s). 
The thick line denotes the limiting case where po ~ 0.0 and 
p= 1.0. 



^ J=lK/L 

j=ii 



(20) 



with ii and ik/l denoting respectively the beginning 
and ending of segment i. The order parameter M2 thus 
effectively tracks the differences in densities of individual 
segments, leading to a measure for the inhomogeneity in 
the system. 

Note that M2 is calculated as an average over T^^^ 
time steps, just as was done for the previous order 
parameter Mi. This means that we calculated the 
summation in equation H19|) at each time step, and 
averaged all these results over a simulation period of 



From the graph of M2 with pq — 0.0 and p — 1.0, it is 
interesting to notice that this order parameter performs 
exquisitely well at tracking the transitions between the 
different traffic regimes. Each time a transition point is 
encountered, M2 drops to zero. This specific behavior of 
the order parameter can be understood by looking at the 
tempo-spatial conditions that hold at these transitional 
densities: each time, the system settles in a complete ho- 
mogeneous equilibrium (as can be seen from the lower 
time-space diagrams in Figs. ITU IT^ and . The dis- 
tribution of the vehicles' space gaps (Fig. also has a 
variance of zero at the first transition point, but not at 
the others. We thus conclude that this inhomogeneity 
measure is able to accurately track the four previously 
discussed phase transitions. 
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V. COMPARISON WITH EXISTING 
LITERATURE 

A final word should be said about the behavior of 
the four phases discussed in this paper. Although the 
study is based on a well-known traffic flow model (i.e., 
the VDR-TCA model), it seems that this behavior 
is drastically different from that observed in real-life 
trafhc. The research elucidated in this paper, is therefore 
relevant as it can be compared to phase transitions in 
other types of granular media, described by cellular 
automata. 

Relating our findings to similar considerations in 
literature, we note the following aspects: Grabolus |3] 
gives a thorough discussion of several variants of the 
STCA, including so-called 'fast-to-start' TCA models 
that give rise to forward propagating density waves; the 
study refers to the two-fluid theory [2^ as a possible 
explanation of this phenomenon. Gray and Griffeath 
discuss the origin of non-concavity in the {k,q) fun- 
damental diagram, where they propose an explanation 
in which the rear end of a jam is unstable, whereas its 
front is stable and growing. Leveque 9] incorporates a 
non-concave fundamental diagram itself in a macroscopic 
model, thereby resembling night driving. Nishinari et al 
[lO| consider the tempo-spatial organization of ants using 
an ant trail model based on pheromones. In their model, 
the average speed of the ants varies non-monotonically 
with their density, leading to an inflection point in 
the {k,q) fundamental diagram. Zhang questions 
the property of anisotropy in multi-lane traffic, leading 
to strikingly similar non-concave {k,q) fundamental 
diagrams. 

Finally, Awazu investigated the flow of various com- 
plex particles using a simple meta-model based on Wol- 
fram's CA-rule 184, but extended with rules that govern 
the change in speed of individual particles. As a result, 
he classified three types of fundamental diagrams, namely 
two phases (2P), three phases (3P), and four phases (4P) 
type systems. Of these types, the 4P-type complex gran- 
ular particle systems closely resemble the regimes discov- 
ered in this paper. Awazu discusses the types of regimes 
and the transitions between them, using the terminol- 
ogy of 'dilute slugs' for the slow moving jams in the 
DCT-phase (see section irVE2|) . which he calls the 'dilute 
jam-flow state'. Analogously, there are 'advancing slugs' 
in the 'advancing jam-flow state' (related to our DAT- 
phase in section IIV E 3|) and the fourth phase is called 
the 'hard jam-flow state' (related to our HCT-phase in 
section ITVE 411 . 

As stated before, the behavior here is different from 
that of real-life traffic flows. Awazu expects that these 
transitions appear if there are more complex interactions 
between the particles. In the 4P-type systems, attractive 
forces between close neighbouring particles and an effec- 
tive resistance on moving neighbouring particles, lead to 



the realization of the previously discussed flow phases. 

VI. SUMMARY 

In this paper, we flrst showed the behavioral char- 
acteristics resulting from the velocity dependent 
randomization traffic cellular automaton (VDR-TCA) 
model, operating under normal conditions (i.e., pq ^ p). 
Then we investigated the more complex system dynamics 
that arise from this traffic flow model in the exceptional 
cases when po ^ P- This behavior was quantitavely 
compared against the VDR-TCA's normal operation, us- 
ing classical fundamental diagrams and histograms that 
show the distribution of the speeds, space, and time gaps. 

Our main investigations were primarily directed at the 
limiting case where pq = 0.0 and p = 1.0. We dis- 
covered the emergence of four different traffic regimes. 
These regimes were individually studied using diagrams 
that show the evolution of their tempo-spatial behavioral 
characteristics. This resulted in the following classifica- 
tion: free-flowing traffic (FFT), dilutely congested traf- 
fic (DCT), densely advancing traffic (DAT), and heavily 
congested traffic (HOT). Our main conclusions here are: 

• all four phases share the common property that 
moving vehicles can never increase their speed once 
the system has settled into an equilibrium, 

• the DAT regime experiences forward propagating 
density waves, corresponding to a non-concave re- 
gion in the system's flow-density relation. 

In order to track the phase transitions between these 
traffic regimes, we looked at two order parameters: the 
first one Mi was based on correlations between neigh- 
boring cells, whereas the second one M2 was based on a 
comparison of the difference between the system's global 
and local densities. We concluded that, as opposed to 
Ml , this latter order parameter M2 performs exemplary 
at tracking the transitions between the different traffic 
regimes. 

Comparing our results with those in existing literature, 
we conclude that the work of Awazu 0], dealing with a 
cellular automaton model of various complex particles, 
gives the closest resemblance to our four phases. 
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